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Abstract 

The existence of "dispersion-managed solitons", i.e., stable pulsating 
solitary-wave solutions to the nonlinear Schrodinger equation with pe- 
riodically modulated and sign-variable dispersion is now well known in 
nonlinear optics. Our purpose here is to investigate whether similar struc- 
tures exist for other well-known nonlinear wave models. Hence, here we 
consider as a basic model the variable-coefficient Korteweg-de Vries equa- 
tion; this has the form of a Korteweg-de Vries equation with a periodically 
varying third-order dispersion coefficient, that can take both positive and 
negative values. More generally, this model may be extended to include 
fifth-order dispersion. Such models may describe, for instance, periodi- 
cally modulated waveguides for long gravity-capillary waves. We develop 
an analytical approximation for solitary waves in the weakly nonlinear 
case, from which it is possible to obtain a reduction to a relatively simple 
integral equation, which is readily solved numerically. Then, we describe 
some systematic direct simulations of the full equation, which use the soli- 
ton shape produced by the integral equation as an initial condition. These 
simulations reveal regions of stable and unstable pulsating solitary waves 
in the corresponding parametric space. Finally, we consider the effects of 
fifth-order dispersion. 
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The existence of "dispersion- managed solitons", i.e., stable pulsating 
solitary-wave solutions to the nonlinear Schrodinger equation with pe- 
riodically modulated dispersion, its mean value being approximately 
equal to zero, is well known in nonlinear optics. However the pe- 
riodic modulation of other nonlinear integrable wave equations has 
received little attention. Here we investigate the Korteweg de Vries 
equation with a periodic dispersion coefficient whose average value is 
small. The model applies to long gravity-capillary waves and natural 
waveguides for internal waves in a stratified ocean. Using analytical 
and numerical techniques, we show that, for harmonic and piecewise- 
constant profiles of the dispersion modulation, stable pulsating soli- 
tons with small amplitudes exist. They become unstable with the 
increase of the initial amplitude or mean dispersion. 



1 Introduction and formulation of the model 

Classical models which give rise to solitons, such as the Korteweg-de Vries 
(KdV) or nonlinear Schrodinger (NLS) equations, describe uniform nonlinear 
waveguides. For example, the NLS equation is a basic model to describe light 
propagation in optical fibers and other guiding structures iQ, while a well- 
known application of the KdV equation is the description of internal waves in 
stratified fluids, propagating in waveguides which exist naturally in the ocean, 
or can be created in a laboratory ||]. In both these models, a (bright) soliton 
exists under a certain condition: in the NLS equation, the dispersion must be 
anomalous while in the KdV equation, the soliton must have a definite 
polarity. 

The necessity to improve stability, bit rate, and other operation characteris- 
tics for soliton streams in optical fibers, have recently attracted a great deal of 
attention to the technique of "dispersion management" (DM), i.e., transmission 
of solitons in an optical fiber which consists of periodically alternating sections 
with opposite signs of the dispersion, so that the average dispersion of the long 
communication link is nearly equal to zero |@, ^, The corresponding NLS 
equation with a variable dispersion coefficient is no longer integrable by means 
of the inverse scattering transform, and no exact soliton solution is known for it. 
Nevertheless, very accurate numerical simulations show that, to an extremely 
high accuracy, this equation gives rise to solitons with an approximately Gaus- 
sian (rather than the ordinary sech), pulsating shape [§1, ||, 0. The existence of 
dispersion-managed solitons of this type is also strongly supported by analytical 
results produced by the variational approximation developed for the DM model 
in different forms & as well as by analysis of the model transformed into an 
integral equation K |10[ (note that the variational approximation can also be 
applied directly to the integral equation |jll|). 

A remarkable property of DM solitons is that they may exist, in a stable form, 
even when the average dispersion, which is not necessarily exactly zero, but may 
be much smaller than the local dispersion in the alternating fiber segments, is 
normal ^ (recall that the ordinary NLS solitons cannot exist in the case of 
the normal dispersion). It is necessary to mention that, prior to the appearance 
of the DM models with the above-mentioned piece-wise constant form of the 
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dispersion profile, the NLS equation with a harmonically (sinusoidally) modu- 
lated sign-changing dispersion term had already been studied in detail, with a 
conclusion, based on both the variational approximation and direct simulations, 
that stable pulsating solitons do exist in that smoothly modulated model p^ . 

In fact, the NLS equation with the local dispersion subject to strong peri- 
odic modulation belongs to a class of periodic heterogeneous systems, in which 
stable pulsating solitons with nontrivial properties may be expected (in this 
work, we use the term "soliton" without assuming integrability of the corre- 
sponding model). In nonlinear optics, other systems belonging to this class are 
tandem waveguides for optical solitons supported by quadratic (x*-^-*) nonlinear- 
ity, based on alternation of x^^-* and linear segments split-step fiber links, 
in which linear segments alternate with those dominated by the Kerr (x^^-*) 
nonlinearity ||l5|| , layered bulk media, with the x*'"^^ coefficient varying between 
the layers so that an optical beam propagating across the layers has its power 
oscillating around a critical value leading to the wave collapse [p^ , and alter- 
nate nonlinear waveguides, composed of periodically alternating waveguiding 
and antiwaveguiding segments . A remarkable feature, common to all these 
systems, despite their very different physical nature, is robustness of the propa- 
gation modes, and the absence of any apparent instability, even when this might 
naively be expected. 

The identification of this class of models, essentially based on the NLS equa- 
tion, makes it natural to ask whether other ordinary soliton-generating equa- 
tions, if subjected to a periodic modulation of the dispersion coefficient, can 
give rise to the propagation of pulsating robust solitary wave (provided that 
a corresponding model describes a physically meaningful situation). The first 
candidate to be investigated is the Korteweg-de Vries (KdV) equation. In fact, 
the KdV equation with variable coefficients is a traditional object for the appli- 
cation of the perturbation theory for solitons However, the case of periodic 
modulation, and the consequent possibility of the existence of a quasi-stable 
pulsating soliton in this case has not yet been studied. 

We start with consideration of the propagation of weakly nonlinear long 
waves in a periodically inhomogeneous waveguide, governed by the variable- 
coefficient KdV equation: 

ut -\- c{ex)ux -\- e [r{ex)uu3: -\- s{ex)uxxx] = O(e^), (1) 

where e <^ 1 is a basic small parameter. Equations of this general form are 
commonly used to describe the propagation of solitary waves in inhomogeneous 
media, for instance, surface and internal waves in stratified fluids Provided 
that the local nonlinear coefRcient, r, and phase velocity, c, are nonvanishing 
functions of ex, one can introduce a propagation coordinate x and a temporal 
variable as follows: 

l^dx, e^t^l^. (2) 

Then Eq. (^ reduces, in the same approximation at which it is valid, to a 
simpler form, 

% -I- 6uug -\- D{x)u0gg = 0, (3) 
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where the local dispersion coefficient is now 



D{x) = 6s/ (rc^) . (4) 
Two integral quantities, which are frequently called mass and momentum, 

r + co r + co 

M= u{e)de, P = / u^{e)de, (5) 



are exact dynamical invariants of Eq. (^), i.e., dM/dx — dP/dx — 0. Note that 
P may be called, more accurately, wave-action flux. 

Our objective is to consider the case when D{x) is a periodic function, 
and, in particular, when it periodically changes its sign. For instance, as it is 
shown in the Appendix, Eq. (|^) describes the propagation of gravity-capillary 
waves over a periodically-varying-bottom topography near the critical value of 
the Bond number; a similar derivation can be easily carried out for interfacial- 
capillary waves. It is also a possible model equation for the propagation of 
internal waves and Rossby waves against the background of a slowly varying 
shear flow In this latter case, the wave modes are localized near critical 
layers, where the long-wave speed is within the range of the shear velocity. 
Generally speaking, in all these applications, one should also take into regard 
the fifth-order dispersion; nevertheless, we initially conjecture that the higher- 
order dispersion is negligible. 

It is assumed that the variable dispersion coefficient D{x) in Eq. (||) can be 
represented in the following, quite general, form: 

D^SDo{x/T) + Diix), (6) 

where Dq is a periodic function with period 1 (i.e., the actual period of the first 
term in Eq. is T) and amplitude 1. The mean value of the function Dq is 
exactly equal to zero. Following the analogy with the DM models in nonlinear 
optics, we will then refer to the parameter S as the map strength. The function 
Di then represents the local average dispersion, which may also be subjected to 
a long-range modulation (a situation with the average dispersion slowly varying 
along the propagation distance is known for fiber-optic DM as well Since 
we are interested in the case of a sign-changing local dispersion, it is usually 
assumed that < S/2. In fact, without loss of generality, one may set 

S = T = 1. Indeed, it is straightforward to see that Eq. (||) is mapped into 
itself, but with 5 = T = 1, by a transformation to new variables, 

X^^, ^^77^' Dir)^Doir) + S-'D,iTr). (7) 



(S'T)i/3' ^ T' (5r)i/3' 

In section 2, we develop an analytical approximation for solitary- wave solu- 
tions to Eq. (||) in the weakly nonlinear case, when the model can be reduced 
to a relatively simple integral equation. Numerically found solutions to the in- 
tegral equations are also given in section 2. Results of direct simulations of the 
full model are displayed in section 3, demonstrating the existence of both stable 
and unstable solitons. In section 4, we briefly consider effects of the fifth-order 
dispersion, and the paper is concluded by section 5. The appendix gives a short 
account of the derivation of Eq. (^), including the fifth-order-dispersion term, 
for a particular hydrodynamic problem. 
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2 Dispersion-dominated waves 

Following the approach of Ref. ||, [T^, where the fiber-optic DM model was 
transformed into an integral equation, we consider Eq. (|^) with the modulated 
dispersion in the form of Eq. with S = T = 1. We assume that the initial 
values supplementing Eq. (||) arc small, viz., f{9) = u{9,x = 0) = 0{e) (note 
that, as it follows from Eq. (0), this is equivalent to assuming that / — 0{1) 
and D varies on a fast scale). It is then natural to introduce a slow time scale 
T = ex, redefine the initial values, so that 

fiO) ^ eg{0), (8) 

and look for a solution in the form 

u{e, t) = eu(°' {9, t) + e^u^i) (9,t) + -- - . (9) 

Defining the "accumulated dispersion" 

W(x) = rOoix) dx, (10) 
Jo 

the zeroth-order solution to Eq. (||) can be found in the form 

^(0) ^ l^^p [,/,3^(^-)j T{A{9,t)}} , (11) 

where J- and J-'^^ stand for the Fourier transform with respect to x ^-iid its 
inverse. From the initial condition (H), it follows that A{9,0) — g{9) in the 
solution (^l|). Because the function is periodic with zero mean, the functions 
W(x) and m^^^ given by Eqs. ( p^ and ( |Tl] ) are also periodic. 
At first order in e, Eq. (||) is reduced to 



(12) 



with u(i)(6',T = 0) = 0. Defining 

h{k,x,T)^T{u^^^+6u^'^^uf^+D,{T)uZ} , (13) 
Eq. ( p^ with the zero initial value has a solution 

u^^^ = -T^^ (cxp (ik^W) r exp [~ik^W{a)] h{k, cr, T)da\ . (14) 



To avoid the growth of secularities in the first-order solution (14), we must 
impose a condition 

/.I 

exp{-tk^W{x))h{k,x,r)dx^O, (15) 

which can be written in terms of the Fourier transform, A{k), of the function 
A{9) introduced in Eq. (pi]), as 



Ar - iDok^A + 'iik / A{K)A{k - k) exp [SikK^K - k)W{x)] dK^dx^ 0. 

(16) 
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If we introduce a function 

y(^) = / exp [inWix)] , (17) 

then it follows from Eq. ( [l^ ) that the amplitude A{9, r) satisfies an evolution 
equation: 

Ar + SJ"-^ l^ik J V{3kK{K- k))A{K)A{k- K)dKj + Do{T)Aggg ^0, (18) 

with the initial condition that A{6,0) — g{9). 

The kernel (|l^) of the integro-differential equation (|l8|) can be easily calcu- 
lated in the case of the piecewise-constant dispersion, which is similar to the 
standard DM scheme in the fiber-optic communications, |0, ||, |^, 

D,{x)-[-"'^ l/4<X<3/4, 
^ ' 1+1/2, 3/4<x<5/4, ^ ' 

which is repeated with the period 1. In this case, the result is 

T/(k) = (8/k) sin(K/8) . (20) 

For the smooth sinusoidal modulation of the dispersion, with 

Z?o(x) = (l/2)cos(27rx), (21) 

the kernel can be found in the form 

{n\f V4/ 

n— ^ ' 

In general, V{k) ^ 1 + O(k^) as k — > 0, while for k ^ 1, and assuming that 
Dq (x) is sufficiently smooth, the integral can be evaluated using the stationary- 
phase approximation. The main contribution then comes from a vicinity of the 
point where W'{x) = -C'o(x) — 0- Assume that I?o(x) has two zeroes, xi and 
X2, which are symmetric, such that X2 = 1 ^ Xi and Z?g(xi) = —D'q{x2) < 
and W{xi) ~ ~W{x2)- Then, the stationary-phase method yields 



Hereafter, we focus on steady-state solitary- wave solutions to Eq. (|T^), fol- 
lowing the approach of Ref. |jl^ . To this end, we assume that the average disper- 
sion Di is a positive constant, and look for solutions in the form A — A{9 — St) 
with a constant velocity S. Anticipating that 5 > 0, we set S = X^, so that the 
Fourier transform A{k) must satisfy an equation following from Eq. (|l^), 



{X^ + Dik^)A = 3 / V{3kK{K- k))A{K)A{k- K)dK. (24) 

J — oo 

By means of the definition A{k) = X^B{k), one can combine two free parameters 
Di and into a single one, Di/X? . Then, given the number N of Fourier modes. 
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the following iterative scheme is used to solve the discrete truncated version of 
Eq. (§): 

o N/2 

Blr^ - y ' p E Vi3klil-k))A^U^:% (25) 

^ l=-N/2 

ir = iii^si"' . (26) 

Here, the superscript is the iteration number, while the integer argument k is 
written as the subscript attached to the amplitudes B^"^ and A^'^\ and 

N-l 2 

P{A} ^ (^^) ' (27) 

is the momentum (see Eq. (||)) of the discrete field Aj. As for similar solutions 
of the NLS equation in Ref. the shape of the steady solution converges; 

however, the momentum diverges, and hence the wave must be normalized at 
each step of the iterative method. 

In Figs. |l| and |[ localized solutions obtained by means of the numerically 
implemented iterative scheme (^), ( p^ ) are shown, as Di is varied for the two 
forms ( |l9| ) and ( ^ ) of the dispersion modulation. Note that localized solutions 
with positive extrema appear to be only possible for Di > 0. By symmetry, 
localized solutions with negative extrema and negative velocity are then only 
possible for Di < 0. As can be seen, a significant difference between the waves 
for the two dispersion-modulation form is observed only near Di = 0. For 
larger Di, the mean dispersion dominates, and the waves for both types of the 
modulation approach the ordinary sech^ soliton of the constant-coefficient KdV 
equation. For clarity, only the domain —20 < a; < 20 is shown in these figures; 
in larger domains, the oscillations apparent in Fig. near Di — decay to zero. 



3 Direct numerical results 

A crucially important issue is stability of the stationary solitary waves obtained 
in the previous section by means of the perturbation expansion. This could 
be addressed by extending the perturbation expansion to the next order, but 
this approach proves unwieldy. A more straightforward approach, which is 
undertaken here, is to use solutions of Eq. (^) as initial conditions for the 
underlying equation (^), and then follow the evolution in direct simulations. 
With the increase of e, the nonlinearity gets stronger (see Eq. (^)), leading 
to stabilization or destabilization of the steady dispersion-dominated solutions. 
To solve Eq. (|^) numerically, we used a standard method combining pseudo- 
spectral techniques in 9 and fourth-order Runge-Kutta integration in x- 

In this section, we consider only the case when the average dispersion Di is 
constant. Further, results presented here are obtained for a slightly smoothed 
version of the piecewise-constant dispersion modulation (p^). Results for other 
forms of the dispersion modulation, including the sinusoidal form (pl|), are very 
similar. Quantitative results, such as the exact position of the stability boundary 
in Fig. 0, see below, may be affected by the choice of the modulation form. 
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Figure 1: Localized solutions of Eq. 
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obtained for the piecewise-constant 
1. In (a) the shape of the pulse, and in 
(b) the variation of the momentum P are shown as D\ is varied. 
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Figure 3: The stability boundary for soliton solutions to Eq. (H) for the piece- 
wise dispersion modulation (|l^). Stable and unstable solutions are denoted by 
asterisks and triangles, respectively. 



10 



The transformation demonstrates that, given a particular form of the 
dispersion modulation, the ratio Di/S and momentum P uniquely define a 
dispersion-dominated solution. Thus, when investigating the stability of dispersion- 
dominated solutions of Eq. (|^), we only have two free parameters, Di/S and 
P. An appropriate value of e used in the previous section is then e = P^^^. 
For the modulation form (|l^, numerically found stability boundaries of the 
dispersion-dominated solutions are shown in figure ^. 

In classifying the solutions as stable or unstable, we have used a simple 
criterion that, for unstable solutions, in some region the Poincare map of the 
solution's amplitude (maximum value) will undergo sustained mean algebraic 
decay in the variable x- Since we solve Eq. (^) with periodic boundary condi- 
tions in 9, at a late asymptotic stage of the evolution, the decay of the soliton 
will be followed by establishment of a nearly uniform state, in which the maxi- 
mum and minimum values of the field are approximately equal and constant in 
X- Conversely, for stable soliton solutions, we expect that, asymptotically, the 
absolute mean of the Poincare map of the maximum value of the field will be 
much larger than its minimum- value counterpart. Strictly speaking, the location 
of the stability boundary depends also on the computational domain size L and 
grid size A9. However, we have taken care to keep L sufficiently large and A9 
sufficiently small to minimize their effect. The boundary also depends on the fi- 
nal propagation distance, Xf- All the solutions shown in figure ^ used Xf = 400, 
which was found to be sufficiently large to make it possible to conclude whether 
a stable soliton was established, or the solution underwent decay. 

Figure ^ clearly shows that there are two distinct regions of stable solitary- 
wave solutions. The stability boundary corresponds to some curve Pt = P{Di/S), 
with the minimum value of Pb occurring at a critical value Rc of Di/S, where 
for the piecewise-constant modulation Rc « 0.25. Then, the stable region in 
< Di/ S < Rc will be referred to as a "dispersion-dominated" one, while the 
stable region in Di/S > Rc will be referred to as a ordinary region. The numeri- 
cal results suggest that the dispersion-dominated region decreases exponentially 
with the increase of Di/S, while the ordinary region increases algebraically with 
the increase of Di/S. In the dispersion-dominated region for Di = 0, it is ap- 
parent that localized solutions of Eq. (||) are only possible up to a finite value of 
the momentum. For the piecewise-constant dispersion modulation, this critical 
value was found to be P w 0.95. Thus, for the variable-dispersion KdV equa- 
tion both nonlinearity and mean dispersion act to destabilize the dispersion- 
dominated solitons. 

The growth of the ordinary stability region with the increase of Di/S is 
easily understood. In this region, the solitary-wave solutions exist due to balance 
between nonlinearity and mean dispersion, the effect of variable dispersion is 
then simply to modulate the form of the wave. In this region, solutions can be 
investigated using the perturbation theory for the ordinary sech^ solitons jlSt , 
or using methods similar to the guiding-center technique, which were developed 
for solitons in optical fibers |Q . 

In Figs . ^, |5| and ^ examples of solutions corresponding to particular points in 
Fig. H are shown. These three soliton solutions are, respectively, stable, unsta- 
ble, and stable again. The dispersion-dominated solution in Fig. |^ demonstrates 
that some momentum is lost from the localized wave to higher-mode oscillatory 
waves. The subsequent interaction between the oscillatory waves and the soli- 
tary wave causes relatively large variability in the Poincare map of the maximum 
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value of the field; however, both the maximum and minimum values eventually 
set down to mean values not significantly different from their initial values. The 
decay of a localized wave is apparent in Fig. ^, where the mean of both the 
maximum and minimum values of the field can be seen to decrease. Eventu- 
ally, at large enough propagation distance they would be expected to become 
approximately equal in the absolute values and opposite in the sign. At this 
stage, the mass and momentum contained initially in the localized wave would 
be completely transferred to the oscillatory waves, as is sugg ested by Fig. ^h). 

For large values of the average dispersion, as is shown in Fig. ^, which is an 
example of a soliton belonging to the ordinary region, it is apparent that the 
solitary- wave solutions are now very similar to the sech^ solitons of the constant- 
dispersion KdV equation. Here, only a small amount of the momentum from the 
initial wave is lost to the oscillatory waves. Because of the reduced amplitude 
of the oscillatory wave, the oscillations in the Poincare map of the maximum 
field value are now much slower than in Fig. ^. 



4 Effects of fifth-order dispersion 

As was mentioned above, the fifth-order dispersion should, generally speaking, 
be added to the model of the DM type with a small average dispersion (similarly, 
third-order dispersion should be included in the optical DM models [^). The 
main effect on the solitons will be to cause them to disperse. In particular, one 
may expect that, for small values of the average dispersion, the addition of the 
fifth-order dispersion to Eq. (|^) will eventually cause stable solutions to become 
unstable. It is expected that this effect is most important when the third-order 
dispersion undergoes a change in sign, i.e. for Di < S/2 in Fig. ^. For larger 
values of Di, it is expected that the main effect of the fifth-order dispersion 
will be to enhance the generation of oscillatory waves j22|. Therefore, here 
we briefly consider numerically the effect of fifth-order dispersion on a stable 
dispersion-dominated solution. 

With the addition of fifth-order dispersion, Eq. (||) becomes 

+ 6uue + D{x)ugee + S'^uggeee = 0, (28) 

where, as shown for the case of gravity-capillary waves in the Appendix, one can 
assume that the coefficient 6^ is constant, in contrast to the strongly modulated 
third-order-dispersion coefficient. In Fig. the effect of varying this parameter 
on the stable solution from Fig. ^ is shown. These numerical solutions are 
obtained using a simple extension of the pseudo-spectral Runge-Kutta method 
described in the previous section. First, the case of 5^ = 10^** was investigated, 
but it is not shown in Fig. ^ as it turns out to be practically identical to 
5^ — Q. Thus, for sufficiently small , stable solitary waves are established as 
well as in the case 5^ = 0. However, as 5"^ is increased, the solitons disperse 
indeed, leaving oscillatory-wave radiation. Other simulations in the dispersion- 
dominated region confirm that the critical (maximum) value of (5^, at which the 
solitary waves are able to form, decreases as both Di/S and P increase, so that 
the critical value is (5^ = on the stability boundary in Fig. B[ 
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Figure 4: A stable dispersion-dominated solution, corresponding to the point 
P = 0.25, Di/S = 0.05 from Fig. (a) The Poincare map of the maximum 
and minimum values of the field, i.e., the maximum and minimum values shown 
at integer values of the propagation distance x- These plots clearly show that 
stable nearly stationary solitons have been found, (b) The form of the soliton 
at X = (displaced by 0.2) and 400. 



13 



(a) 
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Figure 6: The same as in Fig. |j for a stable soliton solution with P = 0.25 and 
Di = 5 in the "ordinary" region. 




Figure 7: The Poincare map of the maximum values of the field, as obtained 
from the numerical solution of Eq. (|8|) for S'^ = 0, 10"'^, 5 x 10"'^, and 10^^. 
Increasing S'^ corresponds to quicker decay of the amplitude. This shows the 
effect of the fifth-order dispersion on the otherwise stable dispersion-dominated 
solution from Fig. El 
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5 Conclusions 



In this work we have introduced a KdV model with a periodicaUy varying disper- 
sion coefficient that takes both positive and negative values. Using a perturba- 
tive expansion based on integral equations and direct numerical simulations, we 
have shown that for the case when the dispersion undergoes a periodic change 
in sign, stable solitary-wave solutions are possible in a region defined by the av- 
erage dispersion and the initial momentum of the system. Further, it has been 
shown that sufficiently weak fifth-order dispersion, which, generally speaking, 
should be taken into regard when the average third-order dispersion is very 
small, does not destroy these solitary waves. 

It should be noted that only one analytical approach (based on the integral 
equation) to investigate the possibility of solitary-wave solutions has been pre- 
sented in detail in this work. We have also tried two additional approaches. 
The first of these is the balance-equation method. Equation (^) can be derived 
from the Hamiltonian 



which, unlike the mass and momentum given by Eq. (g), is not conserved when 
the dispersion coefficient is variable. Instead, it evolves with x according to an 



The balance-equation approach then demands a choice of a slowly varying ansatz 
for the form of the solitary wave. Using the conservation equations for the mass 
and momentum, and the evolution equation ( ^0|) for the Hamiltonian, evolution 
equations for parameters of the ansatz can be derived. The second approach is a 
numerical averaging procedure, involving successive Poincare maps from direct 
numerical simulations of Eq. (H) in an attempt to generate a stable soliton 
solution by iterations. However, neither of these methods have so far proved 
successful for the variable-dispersion KdV equation. 

The reason for solitary-wave solutions being less forthcoming from the variable- 
dispersion KdV equation than for the variable-dispersion NLS equation may 
have to do with the zero-dispersion limit of these two equations. For the 
KdV equation, it is a singular limit with wavebreaking occurring in finite time, 
whereas for the NLS equation in the zero-dispersion limit, only the expected 
phase modulation takes place. Thus, when nonlinearity dominates over the 
variable dispersion in the KdV equation, as in the limit Di/S — and P ^ 1 
in Fig. 1^ above, one would not expect that stable solutions are possible, as 
the evolution generated by the nonlinearity is much faster than that under the 
action of the variable dispersion. Nonetheless, combining the results presented 
here with the well-known ones for the dispersion-managed solitons in the NLS 
models of long fiber-optic links consisting of alternating segments with anoma- 
lous and normal dispersion, we conjecture that stable pulsating solitary waves 
may be possible in a vast class of nonlinear wave equations subject to strong 
periodic modulation of the dispersion coefficient. 




(29) 



immediate corollary of (^) , 




(30) 
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A Zero dispersion limits of long, weakly- nonlinear 
waves: gravity-capillary waves 

A fundamental requirement for the application of Eq. ([!]) is that the linear 
dispersion relationship between the velocity c and the wavenumber k in the 
long wave limit can be written as 

c = co-ak^ + 0{k-^), (31) 

where cq is non-zero, and where for some special parameter combination, the co- 
efhcient a = 0. In geophysical hydrodynamics, the modulation of the coefficients 
of (|l|) can then be induced by a slowly-varying environment, e.g., flow past pe- 
riodically varying topography. We briefly describe here an example where this 
may occur, namely for the flow of a long, weakly-nonlinear, gravity-capillary 
waves over periodically-varying topography when the Bond number is close to 
1/3. A similar derivation holds for interfacial-capillary waves. 

Consider then the propagation of gravity-capillary waves in water of depth 
h with a horizontal length scale L = h/S and amplitude a = eh, such that 
(5, e <C 1. Let g be the gravity acceleration, p the density of the water, and a 
the coefficient of surface tension. The linear dispersion relationship is 

c = co(l - a{hkf + (3{hkf + 0{k^)), (32) 

where 

o.I(I-b), ,.i(|-|-„^) ,33, 

and B = a j pgh? is the Bond number. Therefore, near B — l/i the leading- 
order dispersion coefficient is approximately zero. This occurs for water of the 
depth h « 1cm, in which case one should generally also include diffusive effects. 
However, we will assume here that these are negligible; also we note that for the 
analogous situation of interfacial-capillary waves, the effect of gravity is reduced 
and correspondingly the water depth may not be quite so small, so that in that 
case it is more plausible to ignore diffusive effects. 

For one-dimensional waves propagating to the right, the evolution of the 
free-surface displacement erj is then governed by the extended KdV equation 

Now assume that, in addition, the depth is a slowly varying function of x, such 
that h = h{ex), and introduce new variables 

X = ex, e = g''^(t- j'^y (35) 
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Further, we put v4 = ry/i* and note that in a slowly varying environment the 
wave-action flux density is conserved ||, Hence, the corrected form of 
Eq. (H) is 

OA 3A dA ^2 u^^^'A ,4„^id^A ^. ^ ^3^6. 
'd^-'Wr^^- -W - ' ^'^^ - Oie^eS^n (36) 

Typically a is 0(1); then, let e = (5^, which leads to the variable-coefficient KdV 
equation governing the wave propagation, 

OA 3A dA ^id^A ^ 

However, B = 1/3 is a singular limit of this equation; in this limit we set a — 
and e = J**, hence Eq. (EgI) becomes 



dx 2hVide ^ 80^ dO^ ' ^ ' 

where now (3 = 1/90. Thus, if the topography and Bond number are such that 
the coefficient a is close to zero, and can change its sign, then a uniformly valid 
model equation including the fifth-order dispersion is 

dA_i^dA_ i,d^ _S^d^^ 

dx 2hVi 80 80^ 90 ^-^^^ 
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